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Abstract. The matrix valued analog of the Euler's hypergeometric differen- 
tial equation was introduced by Tirao in [Tj. This equation arises in the study 
of matrix valued spherical functions and in the theory of matrix valued orthog- 
onal polynomials. The goal of this paper is to extend naturally the number 
of parameters of Tirao's equation in order to get a generalized matrix valued 
hypergeometric equation. We take advantage of the tools and strategies de- 
veloped in [l] to identify the corresponding matrix hypergeometric functions 
nFm- We prove that, if n = m -I- 1, this functions are analytic for |z| < 1 and 
we give a necesary condition for the convergence on the unit circle \z\ = 1. 



1. Introduction 

The importance of the hypergeometric differential equation introduced by Euler 
in 1769 and the hypergeometric function was perceived by famous mathematicians 
as Gauss, Kummer and Riemman. Since then, many others found generahzations 
and apphcations. In fact many of the special functions that appear in mathemati- 
cal physics, engineering and statistics are special cases of hypergeometric functions. 
Every second-order ordinary differential equation with three regular singular points, 
by placing the singularities in at 0, 1, and oo, can be transformed into the hyper- 
geometric differential equation 

2(1 - z)f"{z) + {c-z{a + b+ l))/'(z) - abf{z) = 0, 

where a, b and c are supposed to be complex numbers. If c is not an integer 
we can verify by a direct differentiation of the series that the unique solution of 
this equation, analytic and with value 1 at z = 0, is given by the hypergeometric 
function defined by 

{a)n{h)n 



2F1 = 2^ — — - — , where (w)„ = w{w + 1) . . . {w + n ~ 1). 

n>0 ^^>" 

Recently, Tirao introduced in [T] the matrix valued analog of the hypergeometric 
differential equation 

(1) z(l - z)F"{z) + {C-z{A + B + l))F'(z) - ABF{z) = 0, 

where A,B,C € C^'' and F denotes a function on C with values in C. The 
corresponding matrix analog of the Gauss' hypergeometric function is given by 

n>0 



1991 Mathematics Subject Classification. 33C45, 22E46. 

Key words and phrases. Hypergeometric function , matrix spherical functions, matrix orthogonal 
polynomials, fundamental set of solutions. 



2 



P. ROMAN AND S. SIMONDI 



where ( ) o ^^'■^ identity matrix and 

( V )„+! = + ")"'(^ + n){B + n){ %^)^ , for all n > 0. 

The matrix valued hypergeometric functions have connections with the theory of 
matrix valued spherical functions and with the theory of matrix valued orthogonal 
polynomials ([1|,[3],[1])- In fact in [2] we prove that the matrix representations of ir- 
reducible spherical functions associated to the dual symmetric pairs (SU(2, 1), U(2)) 
and (SU(3), U(2)) are given in terms of matrix hypergeometric functions. 

The goal of this paper is to generalize the matrix valued hypergeometric equation 
^ by extending naturally the number of parameters as it was done in the scalar 
case at the end of the nineteenth century. In this way we shall consider the following 
generalized matrix hypergeometric equation 

(2) .A (.A + B,-l)[.± + B.-l)... (.A + B,„ - l) Fi.) 

zf-^A,) (4+A.)...(4 + A„)fW = 0, 

where n, m G Z>o, Ai, . . . , An, . . . , -B,n e C^*" and F denotes a complex func- 
tion with values in C. Observe that if we set n — 2 and m — 1 we obtain the 
Tirao's hypergeometric equation. 

If no eigenvalue of Si, ... , B„i is in the set {0, —1, —2, . . .} we introduce the 
generalized matrix valued hypergeometric function in the following way 



n^m yBi;... ;B^ ' ^ J ~ j] \ Bi; 



3=0 



where the symbol (^^\^: /g" ^ is defined in ([8]). We prove that this functions are 

solutions of the differential equation ([2]). Furthermore if n = m + 1 and {Vj}j-i is 
a basis of C then the set 

( Bi; ...';b'^^ ' ^ }^._^ ' 

is basis of the space of all solutions of the hypergeometric equation Q analytic at 
z = for \z\ < 1. 

Throughout this paper for any matrix A G cx'' we denote by (t{A) to the set 
of all eigenvalues of A and 

p{A) = max{Re{X) : A G cr(A)}, 5{A) = min{Re{X) : A e (t{A)}. 

For any A G C^'', we shall consider the spectral norm || • || defined by 

= max : A is an eigenvalue of ^^*}. 

Let us assume that (cr(Bi) Uct(B2) U . . . Uct(B™)) n (-No) = and that Bi,...,Bm 
are unitarily equivalent to a diagonalizable matrix with real eigenvalues. If 



<5(Bi) + ... + 5(B„0 > IIAill + ... + IIA 



m+l| 



in Thcorem l2.4l we prove that m+iFm si' -bI^^ ' -^j absolutely convergent on 
the unit circle \z\ — 1. 
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2. On the generalized matrix valued HYPERGEOMETRIC EQUATION 

The aim of this section is to give a generahzation of the matrix valued hy- 
pergeometric equation to the case of an arbitrary number of parameters. Let 
Ai, . . . , An, Bi, . . . , B„i G C^''. We begin by considering the following matrix 
valued differential equation of degree max(n, m). 

Let us denote by A to the differential operator z-^. Then, in terms of such operator, 
the matrix valued differential equation ([3]) is given by 

(4) A(A + Bi - 1)(A + B2 - 1) . . . (A + S„ - l)F(z) 

= z(A + Ai)(A + A2) . . . (A + An)F{z). 

We will seek for solutions which have the following form 

00 

where Fj e C. First at all we observe that 



, oc 00 



dz 

3=0 j=o 



and therefore we have 



j=o j=0 



p+j 

3" 



for all 1 < z < n and < fc < m. If we substitute the previous results in the 
differential equation ([4]) we obtain the following recursion relation for the C-valued 
coefficients Fj of F{z) 

(5) {p + ] + l){Bi +p + j)...{B„,+p + j)F,+i 

^ {Ai+p + j)...{An+p + j)Fj. 

If we set j = — 1 in the recursion relation we obtain 

(6) p{Bi + p - 1)(B2 + p - 1) . . . (S„ + p - l)Fo = 0. 
Therefore we have the following indicial equation 

(7) p'' det(Bi + P-1) det(B2 + p - 1) . . . det(B,„ + p - 1) = 0. 

Let Pi, . . . , fS-j. be the eigenvalues of Bj . Then the roots of the indicial equation are 
given by 

p = 0, 1 - . . . , 1 - 1 - . . . , 1 - /32, . . . , 1 - /3r, . . . , 1 - 

The case p — correspond to the analytic solutions of the differential equation ([3]). 
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First at all we will describe the analytc solutions of the equation ([3]). Let 
be the symbol defined recursively by 



(8) 



( bI; ■::-t ) . ^ = (s- + ■ • • (^i + ^r^M (a„ + j) ( 



Ai; 
Si; 



Definition 2.1. If Ai, . . . , An, Bi, . . . , £ i^^rxr eigenvalue of Bi, . 

is in the set {0, —I, —2, . . .}, then we define the funtion 



(9) 



Theorem 2.2. T/ie function 'b\'- "' .'g" ; is analytic for all z £ C if n < m 

and it is analytic on \z\ < 1 if n — m + 1. On the other hand if n > m + f the 
series is not absolutely convergent unless \z\ — 0. 

Proof. Let || • || be any matrix norm on C^''. We observe that 



(A^; ... ;A„ \ 



E 

3=0 



Ai; ... ;A„ 



< 



oo 

E 

3=0 



■Ar, 



Now we will prove that the power series 



3=0 















3 


3<. 



this way we denote Oj 
Observe that 

Gj + l _ \z\ 



is convergent. In 



Then we shall compute lim,_,oo 



< 



.1 + 1 

\z\ 



Ai: ... ■.A„ 
Bi; ... -Bm 



J + 1 



liBn 



J + 1 

■n — m—1 



1 



1 - 



i 



(Bi+jr'(/ii+i) 

-1 



I 



{An+j)\ 



4^1. ..(1 + ^ 



Provided that n < m, this expression tends to zero as j — > cxd and thus the series 
is not absolutely convergent unless z = 0. Ifn>m+1, the series is convergent for 
all z. On the other side, for n = m + 1 the series converges for |z| < 1. □ 

It is worth noticing that if Fq e C then nFm (^^1' " .'g"^ ',z^ Fq is a solution 
of the differential equation In the next theorem we summarize our results 

which characterize the analytic solutions at z = of the generalized hypergeometric 
equation 

Theorem 2.3. Let us assume that (cr(Bi) U a{B2) U . . . U cr(5„i)) n (-Nq) = 
and let Fq £ C then F{z) = „F,„ 

B-i- ... :S„" ; ^) -^0 is a solution of the differential 
equation ([3]) such that F{Q) — Fq. Conversely any solution F analytic at z = is 
of this form. 
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For any A G C'"^'', we shall consider the spectral norm || • || defined by 
\\A\\ = max is an eigenvalue of ylyl*}. 

Since the spectral norm is a unitarily invariant matrix norm we observe that 
= ||A|| for any A e C^'^ and any unitary matrix X e C'^^''. 

Theorem 2.4. Let Ai,..., A„,,+i, Bi, . . . , B„ e C^'' such that (cr(Si) U a{B2) U 
... U a{Bm)) n (— No) is empty and let Bi, . . . ,Bm be unitarily equivalent to a 
diagonalizable matrix with real eigenvalues. If 

S{Br) + . . . + SiB„,) > + + 

then the matrix hypergeometric function m+iFm (^^Bi- '^"^^ ! zj is absolutely con- 
vergent for \z\ = 1. 

Proof. Observe that since Bi is diagonalizable, there exists Xi e C^'" such that 
XiBiX* = Ai where Aj is a diagonal matrix with real no zero entries. Then we 
have that 

\\B-'\\ = \\X,B-'x:\\ = \\Ai\\=p{B-') ^ 



SiBiY 



and 



|(B, + fc/)-i| 



1 



S{B, + kl) 5{B,) + k 
By hypothesis, there exists a positive number A such that 

m m+1 



and we have that 



m+1 



A;! 



1=1 \h=o / i=i 



Ik ™ 



,,L -| m+1 



j\\)k 



< u|fefcEr=i5(B.)-Er=t' 



11^.11 



i=l 



m+1 

n 



i5{Bi))k 

{\\M)k 



If we set \z\ = 1, then 



lim k^+^ 

k^oo 



1 Mi; 
A;! V 



Br, 



J ... 



< lim k 

k^oo 



n 



7r? + 1 



n 



=1 Kik-iy.kW'^ 

(fc__L)!J+^ 
(Pill)fe 



L \{k-l)\k\\^^\\ 



m+1 



= o.nrm)) nr(ll^.ll)-' = o. 
i=l j=l 
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■,A„ 



Since the series X]fc°=i pTx converge if A > 0, the Dirichlet's criterium of numerical 
series of positive numbers implies the absolute convergence on the unit circle for 
\z\ — 1 of the matrix hypergeometric function m+i^m (^^s/- '^^^ ; □ 

Now we will concern ourselves with the non analytic solutions of the matrix 
equation By taking into account the recursion relation ([5]) we introduce the 
following definition. 

Definition 2.5. Let us assume that Ai, . . . , yl„, Si, ... , Bm are r x r complex 
matrices and suppose that -p ^ {a{Bi) + No) U . . . U ((t(S,„) + No) U N. Then we 
define the function 

Observe that p = implies that „F^j ^ "' .'^^ ' ~ ^Fm s^'. 

If /3 e <^{Bi) U • • • U (j{Bm) and (3^1 then p = 1 — /3 is a nonzero solution of the 
indicial equation 

p'' det{Bi +P-1) dct(B2 + p - 1) . . . det(B„, + p - 1) = 0. 

Thus the kernel of the matrix p(i?i + p — 1) . . . {Bm + P ^ 1) is nonempty and we 
can take a nonzero vector G ker(p(i3i + p — 1) . . . {Bm + P — !))• Then, because 
of the way it was constructed, the C-valued function 

F{z) = z^~p^F:^p[j^-:::^:;z)F,, 

is a solution of the matrix generalized hypergeometric equation ^ . We ressume 
this fact in the following theorem. 

Theorem 2.6. Let f3 e cr(Bi) U • • • U cr(S„) and (3^1. If f3 = 1 - p and 

we assume that (3 ^ {a{Bi) + N) U . . . U {a{Bm) + N) U N and Fp is a vector in 
ker(p(Bi + p- 1) . . . {Bm + p - 1)) then 

r(Z)-Z ri^m \ Bt....:B^i^I-^P-/.^p | \Bt+v:...:B„,+vI^^0^ 

is a solution of the generalized hypergeometric equation ([3]) . 

Remark 2.7. In the case that n = m + 1, if the matrices Bi, i = 1, . . . , m, are all 
diagonalizable, and the set of eigenvalues {f3f} is such that (3j ^ if {i, j) ^ (fc, h), 
the set of all analytic solutions of ^ given in Theorem 12.31 and the set of all non- 
analytic solutions given in Theorem 12.61 forms a fundamental set of solutions of the 
hypergeometric equation ([3]). 

3. An example from the representation theory 

The importance of the study of the matrix valued hypergeometric equation, 
and its solutions, become evident at the light of its connections with the theory 
of spherical functions on any iC-type on a Lie group and with the matrix valued 
orthogonal polynomials (See [2l[3l|4j). However the hypergeometric equation that 
appear in [2] is slightly different from equation ([1]) and has the form 

(11) u{l - u)^ + {C- uU)^ - VF{u) = 0, 

ait^ an 
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for r X r complex matrices C, U and V and a C- valued function F on C . This 
differential equation cannot always be reduced to a hypergeometric equation 

(12) z(l - z)F"{z) + {C~z{A + B + l))F'{z) - ABF{z) = 0, 

because we need to be able to find two matrices A and B which are solutions of 
the equations U ~ A -\- B -\- 1, and V — AB, or equivalently we need to solve the 
matrix quadratic equation B^ + {1 — U)B + V = 0. This is not a trivial fact because 
a matrix quadratic equation may not have any solution. For example there is not 
any 2x2 complex matrix X such that 





Now we will concern ourselves in a particular case where it is posible to transform 
the equation (fTT|) into the hypergeometric equation (fT2|) . Let Q{X) = + A(l — 
U) + V. Then det{Q{X)) is a the polynomial in A of degree 2r and therefore there 
are exactly 2r solutions {Ai,...,A2r} of the equation det(Q(A)) — 0. For each 
Ai, we have that the matrix Q{Xi) is singular. Let us assume that we can choose r 
different eigenvalues A^ ^ , . . . , A^^ and r linearly independent vectors pi^, . . . , such 
that pif. £ her Q{Xi^). Let A be the diagonal matrix whose j-th diagonal element 
is given by Xi- and let X be the r x r matrix whose j-th column is the vector pi. . 
Now we observe that 

XA^ + (1 - U)XA - VXA = 0. 

If we multiply the previous equation hy X^^ on the right, we obtain that 

(XAX-^ f + (1 - U)XAX-^ - VXAX-^ = 0, 

and therefore B — XAX^^ is a solution of the matrix quadratic equation + (1 — 
U)B + V = 0. We resume this discussion in the following theorem. 

Theorem 3.1. Let C,U,V G ([^rxr such that we can choose r different eigen- 
values Ai J , . . . , Xi^ and r linearly independent vectors pi^ , • • ■ , Pi^ such that pi^ £ 
kerQ(Aij.). Let A — diag(Ai^.) and let X = [pi^, . . . ,Pi^], i.e. that the j-th column 
of X is Pi - . Then if we let 

B^XAX-^, A = U - XAX-^ -1, 

we have that the differential equation is the following hypergeometric equation 

z{l ~ z)F"{z) + {C - z{A + B + l))F'{z) - ABF{z) = 0. 

Example: Let us consider the differential equation 

. X . ,d'^F(u) , ,dF(u) , , 

(13) u{\-u) +{C-uU)—^-VF{u)^0, 

where the coefHcient matrices are 

i i e 

1=0 i=l i=0 

V = Yi{a + l3 + i-k + l)Eu " - i){i + P-k + l)S,,,+i, 

i=0 1=0 

with a,^ > -1, < A: < /?+ 1 and ^ e N. 
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If we replace a by m G Z>o and by n — 1, this differential equation arises 
from the first few steps in the explicit determination of all matrix valued spherical 
functions associated to the n-dimensional projective space P„(C) = SU(n+I)/U(ri). 
Furthermore, together whith an appropiate choice of a matrix weight W{u), it 
provides examples of families of Jacobi type matrix valued orthogonal polynomials 
(see [3] and [4]). 

The goal is to use the content of Theorem 13.11 to find matrices A and B such 
that U = A + B + 1 and V = AB. First at all we observe that the matrix 
(3(A) = + A(l — U) + V is an upper triangular matrix whose i-th diagonal entry 
is given by 

- \{a + /3 + e + i + 1) + i{a + 13 + i - k + 1). 

Therefore we have that det(Q(A)) = l\''-^Q{X'^-X{a+(3+e+i+l)+i{a+l3+i-k+l)). 
Then we can choose A; as a solution of the cuadratic equation ~ X{a + /3 + £ + 
i + 1) + i{a + (3 + i — k + 1) — and a nonzero Vi £ ker Q{\i) for each Q < i < £. 
For simplicity we shall consider only the case in which A^ ^ Xj for i ^ j- Observe 
that, since Q{X) is an upper triangular matrix, {wi}f=o ^ linearly independent 
set of vectors. Then if we set X = . . . , ve] and A as the diagonal matrix whose 
j-th diagonal element is given by Ajj — Xj, we have that the differential equation 
(fT3| become equal to 

z(l - z)F"(z) + (C - z{A + B + l))F'{z) - ABF{z) = 0, 

where A = U - XKX-^ - 1 and B = XKX-^. 
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